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Abstract 



r^ We present a framework for the analysis of transport processes resulting from the 

^ mean-eddy interaction in a flow. The framework is based on the Transport Induced 

_0 by the Mean-Eddy Interaction (TIME) method presented in a companion paper 

^ [1] . The TIME method estimates the (Lagrangian) transport across stationary (Eu- 

^H ^ lerian) boundaries defined by chosen streamlines of the mean flow. Our framework 

Qh proceeds after first carrying out a sequence of preparatory steps that link the flow 

dynamics to the transport processes. This includes the construction of the so-called 

T—H "instantaneous flux" as the Hovmoller diagram. Transport processes are studied 

^ by linking the signals of the instantaneous flux field to the dynamical variability 

QQ of the flow. This linkage also reveals how the variability of the flow contributes 

to the transport. The spatio-temporal analysis of the flux diagram can be used to 

assess the efficiency of the variability in transport processes. We apply the method 

r^ to the double-gyre ocean circulation model in the situation where the Rossby-wave 

,__! mode dominates the dynamic variability. The spatio-temporal analysis shows that 

T-H the inter-gyre transport is controlled by the circulating eddy vortices in the fast 

^ eastward jet region, whereas the basin-scale Rossby waves have very little impact. 
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1 Introduction 



Analysis of geopliysical flows often employs techniques that decompose the 
velocity field in a manner that will yield a desired insight [2]. A commonly 
used technique is the mean-eddy decomposition 

u(x, t) = u(x) + u'(x, t) , (la) 

Q(x,t) = Q(x) + Q'(x,t), (lb) 

where the field is described by the unsteady eddy activity around a mean state; 
henceforth, {■} and {■} denote the time average (mean) and the residual (un- 
steadiness or eddy), respectively. Here u = (mi, ^2)^ is the velocity field and Q 
represents any property of the fiow such as temperature, chemical or biological 
properties. In the absence of unsteadiness, the kinematic transport occurs only 
along the mean streamlines on which u(x) is everywhere tangent. An effect of 
the unsteadiness is to stir the flow instantaneously and induce mixing across 
the mean streamlines over time. It is also well accepted that instantaneous 
pictures of the unsteady flow themselves do not indicate transport explicitly. 

A variety of Eulerian and Lagrangian methods have been developed to study 
transport observationally, analytically, and numerically. The companion pa- 
per [1] presented a new transport method that is a hybrid combination of 
Lagrangian and Eulerian methods. This paper develops a framework for the 
analysis and diagnosis of transport processes based on this new method. Our 
main focus is on two-dimensional geophysical flows that may be compressible. 
To introduce our method, we begin with a brief discussion of Lagrangian and 
Eulerian methods speciflc to our needs. 

The basis of any Lagrangian method involves the tracking of individual fluid 
particles by solving the initial value problem of J:X = u(x, t). Starting from xq 
at time to? a particle trajectory x(t;xo,to) at time t is given by the temporal 
integral of the local velocity fleld along itself: 



x(t;xo,to) -Xo= / u(x(r;xo,to),r)rfr . (2) 



A very rudimentary description of Lagrangian transport may be obtained from 
a so-called "spaghetti diagram" which is constructed by simply plotting the 
trajectories. Typically this results in a complex tangle of curves from which 
detailed a detailed assessment of Lagrangian transport may prove difficult. In 
recent years the mathematical theory of dynamical systems has provided a 
new point of view and tools for classifying, organizing, and analyzing detailed 



and complex trajectory information by providing a theoretical and computa- 
tional framework for an understanding of the geometric properties of "flow 
structures" . Recent reviews of the dynamical systems approach to transport 
are given in [3,4,5]. 

Nevertheless, there still remain many challenging problems to be tackled by 
Lagrangian methods. Quantifying Lagrangian transport is extremely elaborate 
in general. While techniques based on dynamical systems theory are concep- 
tually ideal for tracking transport of fluid particles, they have not proven as 
useful for studying transport of Q, unless Q is uniform and passive. Moreover, 
Lagrangian methods are not suitable for separating and/or isolating the roles 
played by the mean state u(x) and the unsteadiness u'(x, t) in the trajectory 
x(t;xo,to)- 

In contrast to Lagrangian methods, transport quantities computed with Eule- 
rian methods utilize information taken at pre-selected stationary points. At a 
station x^, the most basic Eulerian transport may be given by the temporal 
integral of the local velocity during a time interval [to,ti]: 



/ u(x^,t)rft=(to-ti)u(x^) . (3) 

Jto 

The resulting transport is associated with u(x), but not u'(x, t) by default 
[compare with (2) for the Lagrangian case]. For transport across a station- 
ary Eulerian curve E = {x-^(p)} where j9 is a parameter along E, the total 
transport during [to,ti] over a spatial segment [pa,Pb] is 



/*' r — x^(p) A u(x^(p), t) dpdt = (ti - to) r -^x^(p) A u(x^(p)) 
Jtn JvA dv JvA dp 



An advantage of the Eulerian methods is the ability to compute the transport 
of Q as well, by replacing u(x, t) with Q(x, t)u(x, t). Overtime, the Eulerian 
methods gives Q{'k, t)u(x, t) = (5(x)u(x)+(5'(x, t)u'(x, t). Hence, the Eulerian 
methods account for the statistical contribution to transport at the second- 
order. 

Our new transport method [1] has the unique ability to identify the effects 
of the mean-eddy interaction in a way that neither Lagrangian nor Eulerian 
methods have accomplished. This advantage that comes from blending of the 
Lagrangian and Eulerian approaches. The method uses information on a sta- 
tionary (Eulerian) boundary curve C to estimate (Lagrangian) transport of 
both fluid particles and Q across C without requiring tracking of individual 
fluid particles. The method estimates the transport by integrating the instan- 
taneous effects of the unsteady ffux while taking the particle advection of the 



mean flow into account. We refer to our metliod tliat quantifies the Transport 
Induced by the Mean-Eddy interaction as "TIME." 

By construction, the TIME method offers a framework for a detailed analysis 
of the spatio-temporal structure of transport processes. The goal of this paper 
is to present this framework through a study of the inter-gyre transport pro- 
cesses in a wind driven, three-layer quasi-geostrophic ocean model (Figure 1). 
Due to its relevance to the mid-latitude ocean circulation, the dynamics of 
wind-driven double-gyre ocean models have been actively studied from vari- 
ous points of view over the last few decades (for review, see [6] and references 
therein). [Fig-1] 

Although the details of transport are highly dependent on the dynamics of the 
flows, there are five common preparatory steps for the analysis using the TIME 
method. The initial two steps of the five concern obtaining an understanding 
of the fiow dynamics based on the mean-eddy decomposition (1). The first 
step is to examine the global fiow structure given by the mean fiow u(x), 
which we call the reference state (Figure la for the wind-driven ocean; see 
Section 2.1 for the details). The second step is to understand the nature of 
the unsteady eddy activity in u'(x, t). Unsteady eddy activity is also referred 
to as the variability. In geophysical fiows, variability is often associated with 
the temporal evolution of the spatially coherent structures (Figure lb; see 
Section 2.2 for the details). It should be clear that these coherent structures 
are defined in the instantaneous Eulerian field and different from the so-called 
"Lagrangian coherent structures" [7]. 

With understanding and insights of the mean and the variability at hand, the 
third step is to compute the the instantaneous fiux that stirs the fiow. The in- 
stantaneous fiux is expressed naturally in terms of a "mean-eddy interaction" . 
(Figure Ic for the wind-driven ocean; see Section 2.3). 

The fourth step is to select the Eulerian boundary, C, of interest based on the 
mean fiow structure. Any mean streamline with reasonable length can be a 
potential C. The actual choice of C should be left up to the specific geophysical 
interests . In Figure la we show our choice of C for the intergyre transport in 
the wind-driven ocean; which we discuss further in Section 2.4). 

The last step is to extract the information of the instantaneous fiux on C. The 
signals of the variability in the mean-eddy interaction are conveniently repre- 
sented by a space-time diagram (i.e., the Hovmoller diagram), which we call 
the flux diagram to emphasize its role in the TIME method (see Section 2.5). 

An important outcome of the five preparatory steps is that they link the fiow 
dynamics to the transport processes, and vice versa, in terms of the mean- 
eddy interaction. This link, and the spatio-temporal integration employed by 
the TIME method, comprise the foundation for the graphical approach to the 



study of the transport processes (see Section 3). In the double-gyre apphcation, 
the analysis will reveal how/when/where the circulating eddy vortices in a 
localized area over the eastward jet are responsible for the inter-gyre transport, 
whereas the basin-scale Rossby waves play a very small role (see Section 4). 
We use the same flow field as in [1] that is nearly periodic in time and has a 
heteroclinic connection in the mean. It is worth noting that the TIME method 
does not require either of such conditions (i.e., presence of the time periodicity 
and the heteroclinic connection). 

The outline of the paper is as follows. Section 2 presents the preparatory steps 
using the application to the inter-gyre transport in the double-gyre wind- 
driven ocean circulation. We extend the TIME functions defined by [1] and 
present a graphical approach that facilitates the analysis of transport pro- 
cesses in Section 3. Inter-gyre transport processes in the double-gyre ocean 
are analyzed in Section 4. Section 5 summarizes the results and provides a 
discussion. 



2 Building the links between variability and transport 



In this section we define the five preparatory steps mentioned in the previous 
section in detail and carry them out in the context of the analysis of inter-gyre 
transport in a wind driven double-gyre ocean circulation model. The flow field 
is obtained by numerical simulation of a three-layer quasi-geostrophic model 
with the model parameters chosen to be consistent with the mid-latitude, 
wind-driven ocean circulations [6,8]. As a result of a constant wind-stress 
curl 0.165dyn/cm^ applied at the ocean surface, the basin-scale circulation 
fluctuates almost periodically around the mean state with dominant spectral 
peak at period T ^ ISldays after the initial 30,000-day spin-up from the 
rest. For our analysis time starts after this spin-up. In this study, we analyze 
the inter-gyre transport processes in the top layer. The instantaneous flow 
patterns are given by the streamfunction '?/'(x, t) that is related to the velocity 
u(x,t)byu(x,t) = (-J^, |:)V'(x,t). 



2. 1 Mean flow 



Given n(x), the corresponding streamfunction '?/'(x) is the reference state that 
provides a geometrical structure of the global flow. Figure la shows '?/'(x) 
for the double-gyre circulation in which the axis of the eastward jet divides 
the ocean basin into two gyres. Because the ocean model has gone through 
the pitchfork bifurcation prior to the Hopf bifurcation [6] (in terms of in- 
creasing wind-stress curl) , an asymmetry exists between the cyclonic subpolar 



gyre and the anticyclonic subtropical gyre. At the confluence of the south- 
ward and northward western boundary currents, the cyclonic subpolar gyre 
and anti-cyclonic subtropical gyre form an asymmetric dipole (along xj and 
xn indicated by the diamonds in Figure la). We refer to this region as the 
"dipole region." The eastward jet defined between the center of the subpolar 
vortex and that of the subtropical gyre and carries the mean net transport 
Tp^rp = j-j/igp — -j/jgj = 24028, where ip^^ (< 0) is the minimum ipix.) over the 
subpolar vortex and ip^^ (> 0) is the maximum over the subtropical vortex. 
The asymmetry of the two gyres is measured by the transport difference, [9] 

'^TD = l^spl - IV'trl = 5438. 



2.2 Dynamic variability of the flow in u'(x, t) 



An understanding of the eddy activity (variability) in u'(x, t) is the key for the 
analysis of transport processes using the TIME method. In the double-gyre 
application, the model ocean dynamics is almost periodic in time resulting 
from the Rossby-wave mode [6,8]. In the left panel of Figure 2, the evolu- 
tion of the ocean dynamics is shown by two time series. The net transport 
of the jet NT{t) = |V'sp(^) — V^st(^)|/'^AfT normalized by the mean state net 
transport ipj^T indicates the fluctuation of the jet strength around the mean. 
The transport difference TD(t) = {\ipsp(t)\ — |'?/'st(^) I )/'^td normalized by iprp^ 
measures the fluctuation of the asymmetry between the subpolar and sub- 
tropical gyres around the mean. Here '^sp(^) is the minimum ■?/'(x, t) over the 
subpolar gyre and ipst{t) is the maximum over the subtropical gyre. The two 
time series show that the amplitude of the fluctuation is of the order smaller 
than 0.1 with respect to the mean. This fact is important since the validity 
of the TIME functions requires the fluctuations to be small compared to the 
mean. 

For convenience, we define the period of the fc-th ocean oscillation by T^^' = 
[tgg + {k — 1)T, t3g + kT) starting from t3g when a minimum of NT{t) occurs 
first time after the spin-up; the number in the subscript of t* represents time 
in day from here on. We also define T^''-^' and T'*^'^' as the first and the second 
half of T^^\ respectively. During T^^'^l, NT{t) increases from a minimum to a 
maximum while TD{t) reaches a minimum after about T/4. Conversely during 
T'^'^1, NT{t) decreases while TD{t) reaches a maximum after about 3T/4 from 
the beginning of T'^^l [Fig-2] 

The eddy streamfunction field ■?/''(x, t) associated with u'(x, t) provides an 
instantaneous pattern of the eddy activity. In the double-gyre application, 
'?/''(x, t) shows two types of eddy activity (Figure lb). One is the westward 
propagation of Rossby waves, whose latitudinally elongated structures are es- 
pecially visible in the eastern basin. Typically there are three waves in the en- 



tire basin, although the one in the western basin is distorted around the dipole 
region. The travel time of the wave from the eastern to the western boundary 
is 453(=151x3)days. Each wave has the longitudinal width 333(=1000/3)km 
and travels with the propagation phase speed 2.2(=1000/453)km/day. 

As mentioned briefly in the introduction, for our study of inter-gyre transport 
in the wind-driven ocean, the axis of the eastward jet is chosen to be C. It is 
the separatrix that connects the western boundary to the eastern boundary. 
We will discuss this in more detail in Section 2.4. 

Although the westward propagation of Rossby waves is seen in almost the 
entire ocean basin, the dipole region has the eddy activity that is more en- 
ergetic. The left panels of Figure 3 show the phases of these eddy vortices 
during T'"^' that starts at tg44. At tg44 + T/4(= tgg2) when TD{t) is minimum 
(corresponding to Figure lb), a strongly positive eddy vortex is located near 
xi. This eddy vortex weakens as the center moves from xi along C as shown 
at t944 + T/2(= t*Q2o) when NT{t) is maximum. It intensifies once again as the 
center reaches X2 as shown at t944 + 3T/4(= t^Qgg) when TD(t) is maximum. It 
moves further along C as shown at tg44 for T later (tg44 + T = t^Qgg), until the 
center leaves C near x^ as shown at tgg2 for 5T/4 later (tg44 + 5T/4(= t^^^^)). 
And it continues to make a cyclic rotation around the sub-polar vortex. One 
cyclic rotation takes 2T. The negative eddy vortex located in the west of this 
positive eddy vortex at t944 + T/4 follows the same cyclic motion but T/2 be- 
hind with the opposite phase of NT{t) and TD{t) with respect to the positive 
eddy vortex. [Fig. 3] 

This cyclic rotation of the eddy vortices is far from uniform because the eddy 
vortices tend to pulsate, i.e., when they intensify during T^^'^'^ or T''^'^] (e-g-, 
^944 + T /A or tg44 + 3T/4 in Figure 3), the centers hardly move; when they 
weaken at the end of T^^^'^l or T''^'^^ (e-g-? ^944 or tg44 + T/2), the centers move 
very quickly. As we shall see below, this complexity in the flow dynamics 
strongly influences the inter-gyre transport processes. We emphasize that these 
eddy vortices are not Lagrangian, i.e., particles don't move with them. 

These two types of eddy activity, westward propagation of the Rossby wave 
and cyclonic circulation of the eddy vortices, are synchronized. In particular 
they merge in the eastern part of the dipole region. At t944 + T/4, a negative 
eddy vortex connects to a negative half of the Rossby wave around X3 as shown 
in Figures 1 and 3. This merger occurs every T/2 and is associated with the 
alternating sign of ■?/''(x, t). 



2.3 Flux variability of (f){'Xi,t) 



The instantaneous flux across the mean flow 



(x,t) = u(x) Au'(x,t) (5) 



is explicitly defined in terms of the "mean-eddy interaction" induced by the 
instantaneous spatial interaction of n(x) and u'(x, t). Geometrically, 0(x, t) 
is the signed area of the parallelogram defined by n(x) and u'(x, t) with the 
unit of 0(x, t) being flux per unit time over the length |u(x)|. The amplitude 
of 0(x, t) depends not only |u(x)| and |u'(x, t)|, but also on the angle between 
u(x) and u'(x, t). We refer to the coherent structures associated with 0(x, t) 
as the flux zones. 

In the double-gyre application, the two types of eddy activity in u'(x, t) lead to 
distinct evolution of flux zones in 0(x, t). In the eastern basin where the Rossby 
waves dominate the variability, the direction change of u(x) in both gyres 
breaks the wave structure latitudinally into three flux zones with alternating 
signs (Figure Ic). Near the jet axis, a sequence of the Rossby waves lead to a 
sequence of flux zones with alternating signs; the positive flux zones correspond 
to northward flux, while the negative ones correspond to the southward flux. 
The width and the westward propagation speed of these flux zones are the 
same as those of the Rossby waves in ip'{'K, t). 

In the dipole region where u(x) is stronger and non-uniform while u'(x, t) 
consists of the cyclic rotation of the four circulating eddy vortices, the mean- 
eddy interaction is complex (Figure 3). Most significantly, the eddy vortices 
along the jet lead to the flux zones with the alternating signs that are visible 
particularly where |u(x'^(s))| is large along the mean jet axis. These flux 
zones in 0(x, t) pulsate in sync with the eddy vortices in ip'{'x., t) (compare the 
right panels of Figure 3 with the left panels). The shapes of the flux zones 
are more loosely defined than those of the eddy vortices due to the spatially 
nonlinear interference by u(x). As the circulating eddy vortices intensify when 
TD(t) is minimum (tg44 -|- T/4) or maximum {tg^/^ + 3T/4), three relatively 
well-defined flux zones are visible over [xj,xi], [xi,X2], and [x2,xn] between 
the eddy vortex centers with the middle one stronger than the other two. 
In between consecutive intensifications, the flux zones weaken and quickly 
propagate along C as NT{t) reaches minimum {tg^^) or maximum {tgj^^ + T /2). 
The propagation speed and direction of the three flux zones are the same as 
those of the circulating eddy vortices. 

As the eddy vortex leaves the mean jet axis near at xn, a weak fourth flux 
zone is generated over [xn,X3] that propagate towards the upstream of the 
mean flow. Because the eddy vortices and the Rossby waves merge around X3 
over [xn,xs], the flux zones associated with them also merge there. 



2-4 Boundary curve C and parameterization of the reference trajectory 



The TIME method uses the streamhnes associated with u(x) as the Eulerian 
boundaries C across which the transport is estimated. Along C, the flight- 
time s is a natural choice of the coordinate variable because C = {x'^(s)} is 
obtained by solving for ^x'^(s) = u(x'^(s)) with a choice of initial position 
x'^(O). Particle advection along C in the mean flow is referred to as the refer- 
ence trajectory. Starting from x'"(so) at to; the reference trajectory is uniquely 
parameterized by sq — to and can be written as (s, t) = {sq — to + t, t) using 
the flight-time coordinate. 

Some key locations on C are shown in Figure 1. In '?/'(x), xj is the location 
that is far enough from the hyperbolic stagnation point of C on the western 
boundary point in s so that n(xj) becomes non-negligible to induce the in- 
stantaneous flux; xn and xs are the locations where the (meandering) jet axis 
makes sharp turns. In ■?/''(x, t), xi and X2 are the locations where the centers 
of the eddy vortices pause to intensify; xn is where the eddy centers leave 
C; around X3, two types of variability, the circulating eddy vortices in the 
upstream and the Rossby waves in the downstream, meet on C. Accordingly, 
xj, xi, X2, Xn and X3 are the boundary points of the flux zones along C in 
0'(x, t) as the eddy vortices pause to intensify and meet the Rossby waves in 
0(x, t). The flight-time coordinates of xj, xi, X2, x^, X3, and xs are sj = 110, 
Si = 114.5, S2 = 118.5, sn = 129, S3 = 150, and ss = 174.5 in the unit of days, 
respectively, starting with x'^(O) = (2 x 10~'^^km, 1011.8km) located very close 
to the western boundary. 



2.5 Flux diagram 

The HovmoUer diagram [10,11] of the instantaneous flux in the {s,t) space: 

/i^(s,t)=0(x^(s),t) = u(x^(s)) Au'(x^(s),t) (6) 

is fundamental to the TIME method because it contains the stirring informa- 
tion locally and instantaneously extracted from 0(x, t) along C. We refer to 
it as the fiux diagram. At a given instance t, a continuous segment of s with 
lJiP{s,t) > corresponds to a positive flux zone where the instantaneous flux 
goes from the right to the left across C with respect to the direction of in- 
creasing s. The direction of the flux is reversed for fi'"{s, t) < 0. The reference 
trajectory (s, t) = (sq — to + t, t) is a diagonal line going through (sq, to) with 
the unit slope (see the main panel of Figure 2). 

The nature of the signals in //'" (s, t) is dependent on both the system and the 
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choice of C. These signals can be complex, as we shall observe in the double- 
gyre application (Figure 2). Nonetheless, having systematically examined the 
mean n(x) (in terms of '?/'(x); Section 2.1), dynamic variability u'(x, t) (in 
terms of ■?/''(x, t); Section 2.2), flux variability 0(x, t) (Section 2.3), and the 
geographic location of C in the flow field (Section 2.4), the physical interpre- 
tation of the signals in fi'-''{s,t) is straightforward. Any signals can be traced 
back to certain flux zones in 0(x, t), and hence the mean-eddy interaction 
process between n(x) and u'(x, t). 

By the construction of the flux diagram along C, the propagation speed and 
direction of the signals in fi'^'{s,t) are deflned with respect to the particle 
advection along the reference trajectory of the mean flow. Signals with positive 
slopes correspond to the downstream propagation of the coherent structures in 
■?/''(x, t). In contrast, signals with negative slopes are related to the upstream 
propagation of the coherent structures. If the slope is steeper than 1, then the 
propagation speed is slower than the particle advection along C in ip{'x.). 

In the double-gyre application, fi'-^{s,t) is almost periodic with period T 
(Figure 2), i.e., iJiP{s,t) = ijp{s^t + T), due to the periodic dynamics in 
il)'{^,t). Within one period, the positive and negative phases are almost anti- 
symmetric, i.e., j/"{s,t) ~ —fi'-^{s,t + T/2). For s < sj, j/"{s,t) is very small 
because of near-zero u(x'^(s)) and changes the sign in synchrony with T''^'^] 
for /i^(s,t) > and Tl'^-^l for fi^^s^t) < 0. 

The strongest signals in jjp{s,t) are concentrated over the segment [sj, sn] 
where the mean jet in '?/'(x) is fast and the circulating eddy vortices in '?/''(x, t) 
are energetic. Due to the pulsation of the circulating eddy vortices as they 
propagate along C (Section 2.3), the corresponding flux zones also pulsate 
simultaneously over the three consecutive segments, Sa = [sj, si), 5*6 = [si, S2), 
and Sc = [s2, sn), with the alternating signs of /i'"(s, t). The widths of Sa, Sb, 
and Sc are narrow (6.5, 4, and lO.Sdays, respectively) with respect to the 
period of intensiflcation T/2 (75.5days) during which the flux zones hardly 
move. Thus the slopes of the dominant signals in fi^'{s, t) are steep. At the end 
of T'^'^] and T'^'^^, these flux zones weaken and quickly propagate downstream 
to the next segment along C. A positive flux zone over Sa during T''^'^^ connects 
to Sb during T^^^'^l, and then to Sc during T^'^'^^-'^l in a sequence; conversely a 
negative flux zone over Sa during T''^'^' moves over to Sb during T'^+^'^l, and 
then to Sc during T^^+'^-'^], Over the subsequent segment Sd = [sn, Ss) along 
C, weak signals of the circulating eddy vortices are observed in ijf[s,t). The 
slope is negative over Sd because of the upstream propagation of the circulating 
eddies as their centers leave C around xn (Figure 3). 

The westward propagation of the Rossby waves is seen for s > S3. The slope is 
negative and less than 1 because the Rossby waves propagate against the mean 
flow with the propagation speed faster than the particle advection along C. 
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Magnitude of f^'"{s, t) rapidly decays to zero for s > 500 because of extremely 
small |TI(x'^(s))| and |u'(x'-^(s), t)| there (not shown). The change in the prop- 
agation slope over s G [200day, 240day] is mainly caused by the meander of 
C (Figure 1). Over Se = [ss,sn), the signals of the Rossby waves and those 
of the circulating eddy vortices are mixed because the two types of variability 
merge between xn and xs in ■?/''(x, t) (Sections 2.2 and 2.3). 

Accordingly, the (s, t) space can be divided into sub-domains based on the 
types of the mean-eddy interaction. Table 1 summarizes the main domains for 
the double-gyre application. Two types of the variability are associated with 
the two domains; D'^^ where the signals are associated with the circulating 
eddy vortices in the dipole region and D™ where the signals are associated 
with the Rossby wave in the downstream region of the eastward jet. The 
total domain is D^^^ = D'^^ U D™. Temporally periodic variability leads to 
the temporal decomposition of D^" based on T''^', i.e., D^^^ = UD^''\ Another 
useful definition D is transient D{t) that covers the (s,t) space up to the 
present time t. [Tab.l] 



3 The TIME functions and the graphical approach to the analysis 
of transport processes 



The companion paper [1] introduced the TIME functions. The main focus in 
[1] was the mathematical formulation, the validity of perturbation approx- 
imations, as well as the verification of the method by comparing with the 
Lagrangian method based on the lobe dynamics for the case of C chosen as a 
heteroclinic connection of the reference state. This section refines the TIME 
functions for a much more detailed analysis of the transport processes due 
to the variability in the flow. We note that although the TIME functions as 
developed in [1] are based on a perturbation approach in the sense that the 
fluctuating part of the velocity field is "small" compared to the reference state, 
none of the development up to this point in the paper requires this smallness 
requirement-all that has been required is the decomposition of the velocity 
field into a (steady) reference state and a fluctuation about the reference state. 



3. 1 The accumulation function 

The accumulation function along C is defined as 



m'^(s,t;to : ti)= / ^f {s -t + t,t) dr . (7) 

•Ito 
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The left-hand side of (7) denotes the amount of fluid transport that occurs 
during the accumulation time interval [to, ti] evaluated at (s, t). The right-hand 
side of (7) expresses the transport in terms of the spatio-temporal integration 
of /i'"(s—t+r, r) during to < r < ti. Thus it can be thought that accumulation 
is advected, while it may continue to occur, with the reference trajectory going 
through {s,t). The sign of m'-^{s,t;to : ti) corresponds to the direction of 
transport across C; m'-''{s,t;to : ti) > is from right to left across C; the 
direction of the flux is reversed for m'-''{s, t; to : ti) > 0. 

The basic idea of (7) is that the transport can be estimated just for the time 
period of interest [to,ti]. It can be extended to study "when", "where", and 
how" variability of the flow contributes to the transport. This is done by re- 
stricting the integration in (7) to a specific space-time domain D that contains 
the particular signals of interest (Section 2.5). Formally and practically this 
leads to the slight modification to the accumulation function: 



m'^{s, t;D)= / H{s -t + T,T; D)if{s -t + r, r)dT 
where 



f8a) 



,1 ii(s-t + T,T) eD, 
H{s-t + T,T;D) = { (8b) 

otherwise 

acts as a switch to turn on and off the instantaneous flux depending on whether 
or not {s—t+T, r) is in D at time r. For example using D{t) as D, m^{s, t; D{t)) 
is the transient transport by accounting for the accumulation up to the present 
time t (see Table 1 for the definition of D{t)). 

For the analysis of transport processes, there are properties of m'^'{s, t; D) that 
are useful (see [1] for technical details). One is the invariance property along 
the individual reference trajectory by advection of the accumulation: 



m^{s,t;D)=m^{s + 6,t + 6;D) (9) 

for any 6, but with D fixed. The invariance property implies a conservation law 
of the accumulation by the advection. At any (s, t) along a reference trajectory, 
m*" (s, t; D) is independent of t as long as D is independent of t. 

The other useful property is the (piece-wise) independence property of D. By 
breaking up the domain D into L non-overlapping sub-domains Di, . . . , D^, 
the piece-wise independence property implies that the accumulation function 
can be written as 
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C/„ +. n\ _ Y^^C 



771 is 



t;D)=J2m^'{s,t;D'). (10) 



1=1 



For example, in the double-gyre application, the effect of the circulating eddy 
vortices and that of the Rossby wave propagation can be examined sepa- 
rately, while the overall transport is the sum of the two, i.e., m*"(s,t; D^^^) = 
m^{s,t; D*^^) + m'-''{s,t;D™). Transient transport can be also decomposed 
into the sub-domains as long as t is the same for all, e.g., wP{s,t\ D^^^{t)) = 

A key feature of the TIME method is the graphical approach using the HovmoUer 
diagram of fi'"{s, t). By definition (8), transport processes are described by the 
manner in which the reference trajectory passes through the signals in fi'"{s, t) 
over the specific domain of interest D (see Figure2). Conversely, it discloses 
the dynamical origins of m'-'{s,t; D) by relating the signals in fi'-'\s,t) to dy- 
namic activities in u'(x, t) through 0(x, t) (Section 2). Thus, the graphical 
approach unveils the dynamical processes of the mean-eddy interaction that 
are responsible for transport. 

Moreover, the graphical approach can assess the efficiency of the dynamic 
activity u'(x, t) in contributing to m*"(s, t; D). The accumulation is, in general, 
extremely efficient if a signal in fi'-''{s,t) propagates with the same unit slope 
because it keeps accumulating the same signed instantaneous flux. In other 
words, transport processes are the most efficient if the (Eulerian) eddy activity 
in '?/''(x, t) propagates at the same speed as the (Lagrangian) particle advection 
of the mean flow ip{x.) along C. In contrast, if the eddy activity is temporally 
periodic and propagates upstream of the mean flow, then the net contribution 
to the accumulation adds up to zero. This is because the reference trajectory 
cuts across the signals of /i'"(s,t) whose sign alternates periodically, resulting 
in the cancellation of the accumulation. 



3.2 The displacement function 



The second type of TIME function quantifies the geometry associated with 
transport, up to the leading-order in the "size" of the unsteady part of the 
velocity field. To describe the geometry in the two-dimensional flow, we use an 
orthogonal coordinate set (/,t) for x near C, where the arc-length coordinate 
/ = l'^'{s) of x'^(s) is deflned along C by ^l^{s) = |u(x'^(s))|. By taking 
the normal projection x'^(s) of x onto C, the signed distance coordinate r = 
r*^(/'^(s)) is deflned by r'^(/'^(s)) = (u(x'^(s))/|n(x'^(s))|) A(x-x'^(s)), where 
|u(x'-^(s))| 7^ 0. If r'^(/'^(s)) > 0, then x lies in the left side of C with respect to 
the direction deflned by the direction of increasing s. The side of x is reversed 
for r^(/^(s)) < 0. 
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The displacement distance function is defined as 



r^{f{s)^,D) = f_\%';^l , (11a) 

u x'" s 



wliere 



a^(s, t;D)= f H{s -t + T,T; D)e^{s - t + r : s)iJ,^'{s - t + t, r)rfr (lib) 



is the displacement area per unit s and H is defined in (8b). Compressibility 
of the flow is taken into account by 



e^(s -t + T : s) 



exp/r D^u(x.^{a))da} , (lie) 

[Js-t+T J 



SO that the flux that has occurred at {s — t + t, t) may expand or compress 
as it advects to (s, t) with the reference particle advection. If the mean flow 
is incompressible as in the quasi-geostrophic model used for the double-gyre 
application, then m'" (s, t; D) and a^ {s^ t; D) are the same since e*" (s — t + r : 
s) = 1 holds for any s — t + T and s. The technical details can be found in [1]. 

The advantage of r'-''{l,t;D) is that it gives the geometry of the so-called 
"pseudo-lobes" that represent the spatial coherency of the transport. Pseudo- 
lobes are the chain-like geometry of transport defined by the areas surrounded 
by C" = {{l,r) I r = 0} and the curve {{l,r) \ r = r^{l,t;D)} at a given 
time t. A positive pseudo-lobe corresponds to a coherent area defined by 
{(/,r) I < r < r'^(/, t; D)}, while a negative pseudo-lobe is a coherent area 
defined by {(/,t) | r'~^{l,t;D) < r < 0}. The relation of the pseudo-lobes to 
the Lagrangian lobes is discussed in [1]. 



4 Application to the Inter-gyre Transport in the Double-Gyre Ocean 



Having carried out the preparatory steps (Section 2) and defined the TIME 
functions (Section 3), we now use the TIME method to analyze inter-gyre 
transport processes in the the double-gyre circulation. We focus on the follow- 
ing three aspects of transport processes using the specific types of sub-domain 
D for m'-^{s,t; D) (see Table 1 for the definition of these sub-domains): using 
D{t), we study how the variability give rise to the development of the pseudo- 
lobes; using D^''{t) and /^''^(t), we compare the impact of the circulating eddy 
vortices and that of the Rossby waves to the inter-gyre transport; using D^^^ 
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along with D{t) and D^''\ we analyze the inter-gyre transport processes of 
the Lagrangian lobes and examine the impact of variability during T''^'. We 
also use mP{s,t\D) in place for a'-''{l,t; D), with which we describe the ge- 
ometry r'-'^[l,t;D) associated with the transport, particular in terms of the 
pseudo-lobes. 

Figure 4 shows the transient accumulation m^{s^t\D{t)) in the Hovmoller 

diagram from the past up to the evaluation (present) time t as t progresses 

upward. Due to the nearly periodic variability in ■i/''(x, t), m'-''{s,t;D{t)) is [Fig. 4] 

also nearly periodic in t, i.e., m^ls, t; D{t)) = m'-''{s, t + T; D{t + T)). Because 

of the anti-symmetry in the oscillation, mP{s,t\D{t)) is also anti-symmetric 

with respect to a T/2-shift in t, i.e., m^{s,t; D{t)) = -m^{s,t + T/2; D{t + 

T/2)). For s < sj at any t, mP{s^t]D{t)) is nearly zero because hardly any 

accumulation happens in the upstream direction of Sj due to extremely small 

yu'^(s, t) there. For s > ss, m'^{s, t; D{t)) ^ rrf{s + 5,t + 5; D{t + 5)) along 

the reference trajectory (e.g., diagonal line in Fig 4) with 5 > indicates very 

little accumulation there. It suggests that the Rossby waves contribute very 

little to the inter-gyre transport in the downstream direction of Ss. 

By a comparison of Figure 4 with Figure 2, accumulation processes in m'-^(s, t; D{t)) 
are closely related to the evolution of the flux zones in ijp{s,t) (Section 2.5). 
Below, we follow the development of a positive pseudo-lobe in Figure 4 starting 
from T^'^'^l as t progresses. This pseudo-lobe first gains positive accumulation 
over S'a where a positive flux zone pulsates during Tt^'^^l If there is no flux 
in the downstream of S'a (i.e., ijp{s,t) = for s > si), then this positive 
pseudo-lobe will spread over [sj, si + T/2] at the end of T^'^-^'^ because the ac- 
cumulation is simply advected along the reference trajectory; this corresponds 
to the invariance property (9). However, positive accumulation that occurs 
over S'a is canceled shortly after it is advected into S'b where a strong negative 
flux zone pulsates during T'^'^l Hence, a small, narrow (in terms of s) positive 
pseudo-lobe develops mostly over S'a during T^'^'^\ At the end of T^'^-^'^ when 
the positive flux zone over S'a moves quickly to S'b, the positive pseudo-lobe 
follows it. In Figure 5, m*" (s, t; D{t)) at t = ^944 -f- T/2 is shown by the dashed 
hne and this positive pseudo-lobe is indicated by T/2. [Fig. 5] 

During T^^'^l, the positive pseudo-lobe continues to develop over S'b where the 
positive flux zone pulsates with large amplitude. Because this flux zone is 
stronger than any flux zones, the pseudo-lobe slowly spreads into the down- 
stream region by the advection of the accumulation. At the end of T'^'^1 when 
the positive flux zone quickly moves to S'c, the positive pseudo-lobe once again 
follows it. In Figure 5, mP{s,t;D{t)) is plotted at t = tg44 + T by the solid 
line. The spread of the pseudo-lobe is observed over S^- 

During T^^-^\ the pseudo-lobe continues to develop over S'c where the positive 
flux zone pulsates before disappearing at sn. It also continues to spread into 
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the downstream region because the flux zones are much weaker there. Speed 
of the spread increases in the downstream direction of S3 because of a weaker, 
positive flux zone induced by the Rossby waves over S^- At the same time, the 
spread of the negative pseudo-lobe that started in T'^'^1 gradually pushes this 
positive pseudo-lobe from the upstream region. At the end of T^^-'^\ the main 
part of the positive pseudo-lobe is located over Se (Figure 5). 

During T'^'^^, the positive pseudo-lobe first gains accumulation from a positive 
flux zone induced by the Rossby waves over Se, but loses it quickly because 
the negative flux zone induced by the Rossby waves moves into Sc- In the 
meantime, the positive pseudo-lobe continues to spread into the downstream 
region. Once the accumulation passes ss, it simply advects along the reference 
trajectory. Near the end of T^^-'^\ the positive pseudo-lobe reaches the flnal 
form as the entire pseudo-lobe passes 53. 

Our analysis of m'^(s, t; D{t)) therefore reveals that the majority of inter- gyre 
transport occurs over a very limited segment [sj,ss]- Accumulation processes 
are synchronized with the evolution of the circulating eddy vortices, and are 
fed mostly by the three flux zones over 3T/2 in t while they are mostly the 
same signed over the spatial segments S^ through Sc] in the subsequent T/2 
in t, the Rossby waves help form the flnal shape of the pseudo-lobe with 
the width T/2 in s. After 2T in t, the pseudo-lobe moves by the advection 
along the reference trajectory because there are no further accumulation in the 
downstream direction of 6*3. Thus, it takes almost 2T in t to develop a fully- 
grown pseudo-lobe of the width T/2 in s. While developing, the propagation 
speed of the pseudo-lobe is much slower than the particle advection along 
the reference trajectory: over 2T in t, the pseudo-lobe moves about T/2 in s 
starting from Sg,. 

The analysis also reveals that the Rossby waves have very little impact on 

the inter-gyre transport for two reasons. One reason is that the signals of 

^'"{s,t) in D™ are weak (see Section 2.5). The other reason comes from 

the fact that the upstream propagation of the signal in n'-'[s,t) in D™ is 

nearly periodic in time (see Section 3.1). Figure 6 shows the decomposition 

of the transient inter-gyre transport m!" [s,t;D{t)) into the part induced by 

the circulating eddy vortices mP{s,t\ D'^^{t)) and that induced by the Rossby 

waves m^(s,t;T)™(t)) at t = t*g + kT (end of T^), where m^{s,t; D{t)) 

= m'^{s,t; D^^it)) +m'-^{s,t; D™(t)) by the piece-wise independence property 

(10). By the deflnition of D'"^, m'^{s,t; T)™(t)) = for s < Sg means that the [Fig.6] 

Rossby waves impact the inter-gyre transport only in the downstream of S3 

(dash-dot line). The amount of accumulation, m'-''{s, t; D''"(t)), is small except 

over [s3, ss] where the signals of the circulating eddy vortices and that of the 

Rossby waves are mixed in ijp{s,t) (Section 2.5). 

Lagrangian lobes have proven to be extremely useful and insightful in a va- 
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riety of transport studies and they provide precise amounts of Lagrangian 
transport of fluid [12,4,5,3]. In the companion paper [1], it was shown that 
m^{s,t; D'^^^) provides a good approximation to the amount of transport car- 
ried by individual lobes in the inter-gyre transport (see also [13]). Figure 7 
shows mP{s^t] D^") by the dash line at the end of T^^l It is doubly periodic 
in s and t because of the periodic dynamics in '?/''(x, t). 

Using the TIME method, we analyze the transport processes associated with 
the Lagrangian lobes. The transient transport m'^'{s,t; D{t)) in the past up 
to the present time t is shown in Figure 7 by the solid line. By the piece-wise 
independence (10), the difference m^{s,t; D'^^^) — m'-'{s,t; D{t)) corresponds 
to the transport that will occur in the future of t. A significant difference is 
observed mainly for s < sn where m^{s,t] D(t)) does not include the active 
accumulation over [sj, s^] yet. A slight difference occurs over [sn, sj + T] due 
to the flux zones of the Rossby waves there. For s > sj + T, m'-^{s,t; D^^^) 
and m^{s,t; D{t)) are almost indistinguishable, suggesting that no further 
transport will occur in the future over the segment of C. Once again we 
confirm that the impact of Rossby waves in the downstream region is neg- 
ligible for inter-gyre transport. Note that m'-''{s, t; D'^^) ^ m^{s, t; D'^^{t)) and 
m'^{s, t; £)'■") ^ nf{s, t; D™(t)) hold over there as well (see Figure6). [Fig.7] 

To examine how much transport occurs during one period T^^^ of the unsteady 
eddy activity in ■?/''(x, t). Figure 7 shows m'-^(s, t; D''^') by the dash-dot line. 
Over [sj,sj + T] of the length T in s, rrf{s^t;D^'^^) and m^{s,t; D(t)) are 
almost indistinguishable because both include the active accumulation over 
[sj, Sn]. Accordingly, during T^^\ only one positive pseudo-lobe grows close to 
its final form over [sj+T/2, sj+T] with width T/2 in s, but no negative pseudo- 
lobe can grow to the final form. If the period is shifted by T/2, then only one 
negative pseudo-lobe grows into its final form over the same [sj + T/2, sj +T] 
segment. For s > sj -|- T, mP{s^t] D''^') differs from mP{s^t] D{t)) because of 
the accumulation in m^{s,t; D{t)) that has occurred prior to T'^^l 



5 Concluding remarks 



Building on the transport method developed in the companion paper [1], we 
have formulated a framework for the analysis of the dynamical processes that 
influence transport. The transport method, called the "Transport Induced by 
the Mean-Eddy interaction" (TIME), is a hybrid combination of Lagrangian 
and Eulerian transport approaches. Our analysis proceeds by a step by step 
approach. In particular, the steps are to determine the mean flow structure 
of ii(x), determine the dynamic variability in u'(x, t), construct the instan- 
taneous stirring chart 0(x, t) = u(x) A u'(x, t) induced by the mean-eddy 
interaction, choose an Eulerian boundary C = {x'^(s)}, and compute the flux 



diagram jjf'^s, t) = 0(x*^(s), t). The signals in fi'-^{s, t) are define relative to the 
reference particle advection along C, which is a diagonal line in the Hovmoller 
diagram. 



The fundamental constructions underlying the TIME method involve com- 
puting transport, either as accumulation irf{s,t\D) or displacement area 
a^{l,t;D) (which gives displacement distance r'"{l,t; D)), which rely on the 
spatio-temporal integration of fi^{s,t). This provides the two-way link be- 
tween the variability of the fiow and the actual transport processes. It is a 
unique feature of the TIME method that neither the Eulerian nor Lagrangian 
methods alone can provide. These fundamentals also provides a platform for 
a novel graphical approach to the analysis of transport processes. 



While transport is highly system dependent, there are some common features 
that can hold in general that we can understand from our graphical approach 
for the analysis of transport processes. For example, the accumulation is most 
effective if the signals in fi'^'{s,t) propagate with the same unit slope as the 
reference trajectory, i.e., dynamic variability propagates with the particle ad- 
vection in the mean fiow. However, if the dynamic variability is temporally 
near periodic, and the signal has a negative slope, then the net effect would 
be almost zero. This may happen when the wave propagates upstream in the 
mean fiow, like in the westward Rossby wave propagation in the double-gyre 
application. The role of variability in transport is analytically studied in [14] 
and is based on a related spatio-temporal scale analysis. 



We have applied our framework to the analysis of intergyre transport processes 
in the double-gyre ocean circulation where the Rossby-wave mode dominates 
the dynamic variability with a period T. The spatio-temporal analysis shows 
that the intergyre transport is controlled by a complex rotation of eddy vor- 
tices in the fast eastward jet near the western boundary current. The emer- 
gence of the pseudo-lobes is synchronized with the circulating eddy vortices 
in u'(x, t). Pseudo-lobes having alternating signed area emerge every T/2 to 
transport water across the mean jet axis between the subpolar and subtropical 
gyres, while each pseudo-lobe sepends almost 2T over a very limited segment 
in the upstream dipole region to fully develop. During the development period, 
the pseudo-lobes propagate at a much slower speed than the reference parti- 
cle advection. However, once fully developed, they propagate downstream of 
the mean jet axis at the same speed as the reference particle advection. The 
basin-scale Rossby wave has very little impact on the intergyre transport. 
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1 Definition of tlie sub-domains for the double-gyre application. 
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Domain 


Description 


Definition 




^cv 


domain associated with the circulating eddy 
vortices 


{{S,t) \ S < Str} 




^rw 


domain asscoiated with the Rossby wave 
propagation 


{{S,t) 1 S > Str} 




^all 


entire domain 


2)cv y ^rw 




Dlk] 


domain for the k-th oscillation period 


{(5,t)|tGTW} 




D{t) 


transient (with respect to present time t) 


{{s,t) \T<t} 




D^^\t) 


transient (with respect to present time t) dur- 
ing the fc-th oscillation period in u'(x, t) 


{(s,t) \ t <t and r 


G TW} 



Table 1 

Definition of the sub-domains for the double-gyre application. 
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List of Figures 



Wind-driven double-gyre circulation at wind-stress curl 
0.165dyn/cm^: a) nieaan streamline field ipix.) averaged over 
T with contour interval 2000; b) eddy streamline field ijj'{'x,t) 
at t = tg44 + T/4(= tggg) with contour interval 500; and c) 
instantaneous flux field 0(x, t) at t = tl^^^ + T/A{= t^^^ with 
contour interval 50. (see also Figure 3). In all panels, the 
dashed contours correspond to the negative values in all panels 
the Eulerian boundary C for the inter-gyre transport shown 
by a thick solid line In '?/^(x), xj, xn, and xg are shown by the 
diamonds; In ■?/'/(x, t) and 0(x, t), Xi, X2, and X3 are shown by 
the triangles. The flight-time coordinates of xj, xi, X2, xn, xs 
and X3 are sj=iio, si = 114.5, S2 = 118.5, sn = 129, sg = 174.5 
and S3 = 150 days. 



The flux diagram fi'-'{s,t) of the intergyre transport along 

the mean jet axis C as the Hovmoller diagram with t runs 

vertically. The contour interval is 20 (km^ / day^) with the 

dashed lines representing negative values. On the abscissa, the 

locations of sj=iio, s^ = 129, and sg = 174.5 are indicated by 

the diamonds, while si = 114.5, S2 = 118.5, and S3 = 150 are 

indicated by the triangles. In jjp{s,t), D^'^^ are indicated by 

the two horizontal lines (solid) at t = ^944 and t = tg44(= t^ggs)' 

while the boundary between Z)'^^ and D™ is shown by the 

vertical line (dashed) at s = S3. The two diagonal lines (solid) 

are examples of the reference trajectory, i.e., (sq — to + t,t) 

with (so, to) = (sj,900) and (so,to) = (250,900). The left 

panel shows NT{t) (solid line) and TD(t) (dashed line) vs. t. 

The top panel shows |u(x'-^(s))| vs. s with the same abscissa 
.c 



as IX {s,t). 



The four phases of 'tp'{'K,t) (left) and (pix^t) (right) during T'^^ 
at t^44 when NT{t) is minumum t*g^^ + T/4:{= tl^^) when TD{t) 
is minimum, t^^^ + T/2(= t^o2o) when NT{t) is maximum, 
and tg44 -|- 3T/4(= tioss) when TD{t) is maximum with time 
increasing upward. The contour intervals are the same as in 
Figure lb for ip'{'x,t) and Figure Ic for 0(x, t). The solid line 
with the diamonds and triangles is C: from the upstream, the 
diamons show xj, Xn, and xg; the triangles show Xi, X2, and 
X3 as shown in Figure Ic. 



24 



Transient accumulation mP[s,t]D{t)). The contour interval 
is 50 with the dashed contours for the negative values; the 
diamonds show sj, sn, and ss while the triangles show si, S2, 
and S3 from upstream to downstream. A reference trajectory 
is plotted for (s, t) = (sq — to + t, t) that goes through ss at the 
end of T[7] with (so,to) = {ss,tl^4 + T). 

Pseudo-lobes of m'^{s,t; D{t)) at the end of T^^-^l 

(t = t*g + (A; - 1/2)T; dash line) and at the end of Tl'^-^l 

(tgg + kT] solid line) for A; > 1; the figure is made using 

k = 7. The positive pseudo-lobe that starts developing from 

the begining of Tt^-^l at tl^^{= t^g + (k - 1)T with A; = 7) is 

indicated by: T/2 at t = tl^^ + T/2; T at t = t^^^ + T; 3T/2 at 

^944 + ST/2; and 2T at ^944 -|- 2T. Dark diamonds show sj, sn, 

and ss; dark triangles show Si, S2, and S3; lighter diamonds 

show sj -|- T, Sn + T, and ss + T; and light triangles show 

si +T, S2 + T, and S3 -|- T. 1301 

Pseudo-lobes of m^{s,t; D{t)) for the total transient transport 
(solid line; same as in Figure 5), irf{s,t\D'^'^{t)) by the 
circulating eddy vortices (dash line), and 7iiP{s,t;D™{t)) by 
the circulating eddy vortices (dash-dot line) at t = tgg -|- kT 
(i.e., the end of Tl'^l); the figure is made at t^ggg using k = 7. 
Dark diamonds and dark triangles are the same as in Figure 5. [3T] 

Pseudo-lobes of mP{s,t\ D{t)) for the total transient transport 
(solid line; same as in Figures 5), m'"(s,t; -D''^') for transport 
induced during Tt'^l (dash-dot line), Tn'^{s,t; D^^^) for total 
transport (dash hne) at t = t^Q + kT (i.e., the end of T^'^^); the 
figure is made at t^ggg using k = 7. Diamonds and triangles are 
the same as in Figure 5. 
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Fig. 1. Wind-driven double-gyre circulation at wind-stress curl 0.165dyn/cni^: a) 
meaan streamline field ip{x) averaged over T with contour interval 2000; b) eddy 
streamline field ijj'{x,t) at t = ^944 -|- T/A{= ^^2) with contour interval 500; and c) 
instantaneous flux field <j){x,t) at t = ^944 -|- T/4{= ^932) with contour interval 50. 
(see also Figure 3). In all panels, the dashed contours correspond to the negative 
values in all panels the Eulerian boundary C for the inter-gyre transport shown by 
a thick solid line In ip^x), xj, xn, and xs are shown by the diamonds; In 'ipf{x,t) 
and 0(x, t), xi, X2, and X3 are shown by the triangles. The flight-time coordinates 
of xj, xi, X2, Xn, Xs and X3 are sj=iio, si = 114.5, S2 = 118.5, sn = 129, ss = 174.5 
and S3 = 150 days. 
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Fig. 2. The flux diagram /U (s, t) of the intergyre transport along the mean jet axis 
C as the Hovmoller diagram with t runs vertically. The contour interval is 20 (km^ 
/ day^) with the dashed lines representing negative values. On the abscissa, the 
locations of sj=iio, sn = 129, and ss = 174.5 are indicated by the diamonds, while 
si = 114.5, S2 = 118.5, and S3 = 150 are indicated by the triangles. In i^'^{s,t), D'-'^' 
are indicated by the two horizontal lines (solid) at t = ^944 and t = i944(= ^1095)) 
while the boundary between D^^ and Z?''™ is shown by the vertical line (dashed) at 
s = S3. The two diagonal lines (solid) are examples of the reference trajectory, i.e., 
(sq — to + t,t) with (so, to) = (sj, 900) and (sq, to) = (250,900). The left panel shows 
NT{t) (solid line) and TD{t) (dashed line) vs. t. The top panel shows |u(x'^(s))| 
vs. s with the same abscissa as iJp{s,t). 
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Fig. 3. The four phases of ip'{'K,t) (left) and </)(x, t) (right) during T^'^J at tg^^ when 
NT{t) is minumum tg^^ + T/4:{= ^^2) when TD(t) is minimum, tl^^ + T/2{= tio2o) 
when NT{t) is maximum, and ^944 + 3T/4(= ^1053) when TD{t) is maximum with 
time increasing upward. The contour intervals are the same as in Figure lb for 
■(/''(x, t) and Figure Ic for (j){x,t). The solid line with the diamonds and triangles is 
C: from the upstream, the diamons show xj, xn, and xs; the triangles show xi, X2, 
and X3 as shown in Figure Ic. 
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Fig. 4. Transient accumulation mP {s,t;D{t)). The contour interval is 50 with the 
dashed contours for the negative values; the diamonds show sj, sn, and ss while the 
triangles show si, S2, and S3 from upstream to downstream. A reference trajectory 
is plotted for (s, t) = (sq — Iq + t, t) that goes through ss at the end of T'^J with 
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Fig. 5. Pseudo-lobes of m^{s, t; D{t)) at the end of T^^-'^^ {t = q^ + {k- 1/2)T; dash 
hne) and at the end of T'*^'^] (igg + kT; soUd hne) for A: > 1; the figure is made using 
k = 7. The positive pseudo-lobe that starts developing from the begining of T'-'^-^' 
at i944(= ^36 + {k- l)r with k = 7) is indicated by: T/2 at t = tl^^ + T/2; T at 
t = tg44 + T; 3r/2 at tg44 + 3T/2; and 2T at tg44 + 2r. Dark diamonds show sj, sn, 
and ss; dark triangles show si, S2, and S3; lighter diamonds show sj + T, sn + T, 
and Ss + T; and light triangles show si + T, S2 + T, and S3 + T. 
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Fig. 6. Pseudo-lobes of m {s,t; D{t)) for the total transient transport (solid line; 
same as in Figure 5), m {s,t]D'^^{t)) by the circulating eddy vortices (dash line), 
and mP {s, t; L'''^(t)) by the circulating eddy vortices (dash-dot line) at t = tgg + kT 
(i.e., the end of T^'''); the figure is made at tiggs using k = 7. Dark diamonds and 
dark triangles are the same as in Figure 5. 
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Fig. 7. Pseudo-lobes of m {s,t; D{t)) for the total transient transport (solid line; 
same as in Figures 5), m {s,t;D^^) for transport induced during T^ ' (dash-dot 
line), mP {s,t]D'^^^) for total transport (dash line) at t = tgg + kT (i.e., the end of 
T^^); the figure is made at ^1995 using k = 7. Diamonds and triangles are the same 
as in Figure 5. 
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